Microscopic Derivation of Causal Diffusion Equation using Projection Operator 

Method 



T. Koide 1 

1 Instituto de Fisica, Universidade de Sao Paulo, CP. 66318, 05315-970 Sao Paulo-SP, Brazil 

We derive a coarse-grained equation of motion of a number density by applying the projection 
operator method to a non-relativistic model. The derived equation is an integrodifferential equa- 
tion and contains the memory effect. The equation is consistent with causality and the sum rule 
associated with the number conservation in the low momentum limit, in contrast to usual acausal 
diffusion equations given by using the Fick's law. After employing the Markov approximation, we 
find that the equation has the similar form to the causal diffusion equation. Our result suggests that 
current-current correlations are not necessarily adequate as the definition of diffusion constants. 
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I. INTRODUCTION 

Diffusion is a typical relaxation process and appears 
in various fields of physics: thermal diffusion processes, 
spin diffusion processes, Brownian motions and so on. It 
is empirically known that the dynamics of these processes 
is approximately given by the diffusion equation, 



L»V 2 M(x,t) = 0, 



(1) 



where D is the diffusion constant. This equation is phe- 
nomenologically derived by employing the Fick's law or 
the Fourier's law. 

Although the diffusion equation has broad applica- 
bility, there exist the limits of the validity. First 
of all, the diffusion equation does not obey causality 

BiHaiBHSHliamEEQ, the propa- 
gation speed of information exceeds the speed of light. 
This means that we cannot apply the ordinary diffu- 
sion equation to describe relativistic diffusion processes 
that mi ght be realized in relativistic heavy- ion collisions 
[ToL [Til Il2t fl3| . Second, the diffusion equation breaks 
sum rules associated with conservation laws jl 51]. The 
diffusion equation is the coarse-grained equation that is 
valid only for describing macroscopic motions, and hence 
one may claim that such a coarse-grained dynamics does 
not necessarily satisfy the sum rules. However, as we will 
see later, it is possible to derive a coarse-grained equation 
consistent with a sum rule. 

These deficiencies can be overcome by introdu cing re- 
laxation times p, s n a S 1 SS SES IH IH III 13 

liH ITU ITsl ligj . Then, the diffusion equation is changed 
into the following telegraph equation: 

r^u(x, t) + — u(x, t) - DV 2 u{x, t) = 0, (2) 

where D is the diffusion constant and r is the relaxation 
time. The telegraph equation is reduced into the diffu- 
sion equation in the limit of r — > 0. The propagation 
speed of the equation is defined by V = y/D/r. One 
can easily see that the propagation speed of the diffusion 
equation diverges and hence causality is broken. Thus, 



in the following, we call the diffusion equation given by 
the Fick's law the "acausal" diffusion equation, and call 
the telegraph equation the "causal" diffusion equation. 
Furthermore, as we will see later, the causal diffusion 
equation does not break sum rules. 

As just described, the causal diffusion equation may 
be more appropriate to describe diffusion processes [20j. 
However, the microscopic derivation of the causal diffu- 
sion is still controversial. One of the typical methods to 
derive coarse-grained equations is the p ro j ection o per ator 

methodf pomiM El m M M m, 111 IH Sli m 

I32ll35 . l34l35ll36l |. In this method, the motions associated 
with microscopic time and length scales are projected out 
by introducing projection operators, and we can obtain 
coarse-grained equations expressed in terms of variables 
associated with macroscopic time and length scales. As 
is discussed in Ref. [22j , the coarse-grained equation de- 
rived by using the POM is not causal diffusion equations 
but acausal diffusion equations. However, it should be 
noted that an approximation whose validity is not obvi- 
ous is introduced to obtain the acausal diffusion equation. 
Thus, when we do not use the approximation, it may be 
possible to obtain causal diffusion equations instead of 
acausal ones by using the POM. 

As a matter of fact, the author recently applied the 
POM to derive the coarse-grained equation of the order 
parameter that describes the critical dynamics of the chi- 
ral phase transition [23Ll24| . The derived equation fulfills 
the requirements near the critical temperature and con- 
verges to the equilibrium state consistent with mean-field 
results evaluated in finite temperature field theory. How- 
ever, the equation shows the relaxation exhibiting oscil- 
lation. This behavior is different from that of the time- 
dependent Ginzburg-Landau (TDGL) equation that has 
been assumed as a phenomenological equation of the crit- 
ical dynamics. We can look upon the TDGL equation as 
the acausal diffusion equation for " non-conserved" quan- 
tities because of its overdamping behavior. Then, the 
appearance of the oscillation means that a kind of relax- 
ation time is introduced. This result suggests to us that 
if we apply the POM to derive a coarse-grained equation 
of a "conserved" quantity, for instance, a number den- 
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sity, the coarse-grained dynamics may be accompanied 
by oscillation like the causal diffusion dynamics. 

In this paper, we apply the POM to a non-relativistic 
model and show that the coarse-grained equation of the 
number density has a similar form to the causal diffusion 
equation instead of the acausal one. 

This paper is organized as follows. In Section 2, we 
summarize the POM. In Section 3, we apply the POM 
to a non-relativistic model and derive the coarse-grained 
equation of the number density. The number density is 
a conserved quantity in this model and there exists the 
sum rule associated with the conservation law. We inves- 
tigate the relation between the coarse-grained equation 
and the sum rule in Section 4. The coarse-grained equa- 
tion contains the memory effect, which can be eliminated 
by employing the Markov limit. Then, the equation has a 
similar form to the causal diffusion equation, as is shown 
in Section 5. In the acausal diffusion equation of the 
number density, it is known that the diffusion constant is 
given by the time correlation function of the number den- 
sity. However, the simple relation is changed for causal 
diffusion equations. The reason is discussed in Section 
6. The summary and concluding remarks are given in 
Section 7. 



II. PROJECTION OPERATOR METHOD 

In this section, we giv e a short review of the POM 

The time-evolution of an arbitrary operator follows the 
Heisenberg equation of motion, 



where L is the Liouvillc operator and to is an initial time 
at which we set up an initial state. To carry out the 
coarse-grainings of irrelevant information, we introduce 
the projection operators P and its complementary oper- 
ator Q = 1 — P with the following generic properties: 



P = P, 
PQ = QP = Q. 



(6) 
(7) 



From Eq. one can see that operators are evolved by 
e iL(t-t ) ^ fast obeys the following differential equation: 



jl e i£(t-to) 
dt 



e iL{t-t ) iL 

e ^(t-t ) piL + ^L(t-t )Q iL ^ 



Operating the projection operator Q from the right, we 
obtain 



JL(t-t ), 



o iL(t-t ) 



dt 



PiLQ 



o iL(t-t ) 



QiLQ. (9) 



This equation can be solved for e lL{ - l ^Q, 

e iL(t-t )Q 

ft 

Q e iLQ(t-t ) 



f dTe lL(T - to) PiLQe lLQ(t ' T) .(10) 

J to 





= i[H,0(t)] 


(3) 


/' 

J to 


= iLO(t) 


(4) 


where 


^ 0(t) 


= e lL(t -* 0) O(t ), 


(5) 


V(t,t, 



Substituting Eq. (|10|> into the second term on the r.h.s. 
of Eq. JHJ and operating O(to) from the right, the Heisen- 
berg equation of motion is rewritten as 

j t O(t) = e lL{t - to) PiLOfo) 

+ f dTe lL(t " T) PiLQe lLQ(T - to) iLO(t ) 

J to 

+Qe lLQ{t - to) iLO(t () ). (11) 

This equation is called the time-convolution (TC) equa- 
tion. The first term on the r.h.s. of the equation is 
called the streaming term and corresponds to a collec- 
tive oscillation such as plasma wave, spin wave, and so 
on. The second term is the memory term that causes 
dissipation. The third term is the noise term. We can 
show that the memory term can be expressed by the 
time correlation of the noise. This relation is called 
the fluctuation-dissipation theorem of the second kind 
[ll|2l|2l|^|M|M|3l|3l|3l|M|36|. However, in 
this paper, we simply drop the noise term in the following 
discussion. 

The TC equation l|ll|l is still equivalent to the Heisen- 
berg equation of motion and difficult to solve in general. 
Note that we can reexpress the memory term of the TC 
equation as 



dTe lL{t - T) PiLQV(T, t Q )e lQLoQ{T - to \LO(t Q ), (12) 



1 + J2 l " / d< i 



ti 



tn-l 



with 



dt 2 - ■ ■ / dt n 

n— 1 u *° ^ *° 

xQL?(t n - i )Q£?(i„-i - t Q ) ■ ■ ■ QLf(h - to), 

(13) 

lf(t-t ) = e 4( 3 i »Q( t -*«)L /e - l(3ioQ( *- t ^. (14) 



Here, Lq and L/ are the Liouville operators of the non- 
perturbative Hamiltonian Hq and the interaction Hamil- 
tonian Hi, respectively, 



L O = [H ,O], L I = [H I ,0}. 



(15) 



When we expand V(t, to) up to first order in terms of Lj, 
the TC equation is given by 

4o(< ) = e^-^PiLOito) 



dt 



+ / dTe lL{t - T) PiLQe lQLoQ{T - to) iLO{t ). 

J t 

(16) 



Here, we have already dropped the noise term. Equation 
(|16fl is the starting point in the following calculation. 

Note that the memory term contains the coarse- 
grained time-evolution operator e l Q L oQt^ j n general, it 
is not easy to evaluate the coarse-grained time-evolution 
and hence it is approximately replaced by the ordinary 
time-evolution operator e* Lot p^ . 

f dTe lL{t - T) PiLQe lQLoQT iLO{Q) 
Jo 

w [ dre lL{t ~ T) PiLQe lLaT iLO{0). (17) 
Jo 

See Appendix^for details. This approximation is crucial 
to derive the "acausal" diffusion equation. However, as 
we will see later, it yields several problems. 



III. APPLICATION TO NON-RELATIVISTIC 
MODEL 

In this section, we apply the POM to a non-relativistic 
model. In the view of causality, coarse-grained equations 
of non-relativistic models is not necessarily causal diffu- 
sion equations. However, if we can implement coarse- 
grainings preserving conservation laws, the causal dif- 
fusion equations are more appropriate as coarse-grained 
equations even in the non-relativistic systems, because 
acausal diffusion equations break sum rules. 

We apply the POM to the non-relativistic model with 
the following Hamiltonian: 



H = 
H = 

Hi = 



Hq + Hi, 
d 3 x^ f ( x ) 



2m 



(JL V(x), 



(18) 
(19) 



- / d 3 xd 3 x'|^ t (x)V'(x) W (x-x')V t (x')V'(x'), 

(20) 



where Hq and Hi are the non-perturbative Hamiltonian 
and the interaction Hamiltonian, respectively. The chem- 
ical potential /i is introduced in the non-perturbative 
Hamiltonian. 

The commutation relation of the fermion field ^(x) is 
given by 



[V(x),^(x')] + =^ 33 (x-x0, 
- denotes the anticommutator. 



(21) 



where [ 

There are many possibilities for the definition of the 
pr ojec tion operator |22l l23t 124 El El El MM MM 
|3l|3!|3E!E3 ' ^° d escrH3e the diffusion process of the 
number density, we employ the following Mori projection 
operator |27|. 



PO 



d 3 xd 3 x'(0,Sn(x)) • (<yn(x),,5n(x')) _1 • Jn(x'). 



(22) 





FIG. 1: The ring diagram that contributes the calculation of 
the memory term. The solid line represents the fermion and 
the dotted line denotes the interaction w(x). 



Here, 5n(x) denotes the fluctuations of the number den- 
sity, 

fe(x) = ^(x)V(x) - (^ t (x)V(x)} e9 . (23) 

where (ip* (x)ip(x)) eq is the expectation value in thermal 
equilibrium. The Kubo's canonical correlation is defined 
by 

r?d\ 
(3 



( X , Y ) = / -7 



{ e XH °Xe- XH < 



°Y) C 



(24) 



where 



(O)o 



1 



Ti[e- pH "0] 



(25) 



with Zq = Tr[e~" H °]. When we use this projection op- 
erator, we obtain a linear equation. This is enough for 
our purpose because we are interested in the dynamics 
close to equilibrium. If we want to take nonlinear terms 
into account, we should re-define a projection operator 
including nonlinear operators. See |24l l27j for details. 

Substituting them into the TC equation (|lfi[) and set- 
ting O(0) = Sn(x) and to = 0, the coarse-grained equa- 
tion of the number density is derived. As is discussed 
in Refs. the memory function is approximately 

given by calculating the contribution of the ring diagram 
shown in Fig. ^ Then, we set v(x — x') = Jw(x — x') 
for simplicity. Finally, we have the integrodifferential 
equation of the number density, 



dt 



i(k,t) = - [ dsT(k,t 
Jo 



s)5n(k, s) 



(26) 



Here, the memory function T(k, t) is defined by the 
inverse-Laplace transform of r L (k, s), 

*(k) 



where 



r*(k,a) 

fi(k) 
Xf(k) 



-X s 



X0 (k)+xf(k) 



n(k), 



= 1 - 5/3X0 (k), 



dte 



dt 



Xt(k), 



r°° d 2 

(k) = j o dte~ st — Xt (k), 



(27) 

(28) 
(29) 

(30) 
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with 



Xt(k) 



-i(E p 



p+k 



Here, n k is the Fermi distribution function, 

1 



3 /3(B k -M) 4. 1' 



(31) 



(32) 



with — k 2 /2m. The function x*(k) corresponds to 
the simple contribution of the one-loop fermion diagram 
and is interpreted in the kinematical way as is discussed 
in Refs. [23,|24(. However, the memory function itself 
is given by the involved combination of Xt(k) because 
of the coarse-grained time-evolution operator. If we use 
the approximation discussed in Eq. l|17fl . the memory 
function is more simplified. We will come back to this 
point later. 

The Laplace transform of the diffusion equation i|26|) 



Sn L (k,s) = 



5n(k, 0) 

s + r L (k,s)' 



(33) 



Now, we can investigate whether the number density con- 
verges to the thermal equilibrium state or not. When 
the diffusion equation describes the thermal equilibration 
process, 5n(x, t) vanishes at late time. From the final 
value theorem of the Laplace transformation, <5n(k, 00) 
is given by 



lim Sn(k, t) = lim sSn (k, s). 

t — >oo s — >0 



(34) 



Substituting Eq. i|33|) into the final value theorem, we 
have 



lim Sn(k, t) = 0. 

t— *-oo 



(35) 



Thus, the time-evolution of 5n(k, t) is consistent with 
the fact that the derived diffusion equation describes the 
thermal equilibration process. 

As we have discussed at the last paragraph of the pre- 
ceding section, the projection operator Q contained in 
the memory term is sometimes dropped and the approx- 
imation yields several problems. One is the problem of 
the convergence discussed above. When we apply the 
approximation, from Eq. (|17|l . the memory function is 
given by 



r L (k, s ) 



#(k) 
Xo(k) 



fi(k). 



(36) 



Substituting this expression into the final value theorem, 
one can easily check that <5n(x, t) does not vanish at late 
time, 



lim 5n(k, t) ^ 0. 



(37) 



This means that the derived equation cannot describe the 
thermal equilibration process in employing the approxi- 
mation. Thus, we should not apply the approximation 
in this calculation. Another problem will be discussed in 
Section 6. 



IV. COARSE-GRAININGS AND SUM RULE 

In our Hamiltonian, the number density is a conserved 
quantity and there exists a sum rule associated with the 
conservation law. Then, it is desirable that the coarse- 
grained equation is consistent with the sum rule. In this 
section, we show that our coarse-grained equation is con- 
sistent with the sum rule. 

First of all, we introduce the Fourier transform of the 
correlation function of the number density |l5j , 



<[n(x,t),n(x',t')]>e 



du 



dk 



2vr J (2tt) 3 



C"(k,w)e 



ik(x— x') iLo{t — t ) 



(38) 



The Fourier transform C"(k, lo) is real and an odd func- 
tion of the frequency u> |15|. From the sum rule asso- 
ciated with the number conservation, as is discussed in 
Appendix^] the Laplace- Fourier transform (5n LF (k, z) is 
expanded for large values of z, 

Sn LF (k, z)/F{k) 

= l -C(k) + ^-k 2 (n(0)) eq + O(l/z% (39) 
z z a m 

where F(k) denotes an external field, and 

du' C"{k,u)') 



C(k) 



2tt 



(40) 



We can see that (i) the term proportional to 1/z 2 dis- 
appears and (ii) the coefficient of the term 1/z 3 is given 
by i(n(0)) eq k 2 /m. As is shown in Appendix iBl if the 
dynamics of the number density is approximately given 
by the acausal diffusion equation, one can show that the 
term proportional to 1/z 2 does not disappear |15|. This 
is contradiction to the exact result. 

On the other hand, the integrodifferential equation 
(|26|1 obeys the sum rule and is consistent with Eq. 13911 . 
By setting s = —iz and expanding Eq. l-i-il) for large 
values of z, we obtain 



Sn LF (k,z) 



iSn(k, 0) 

z 

i5n(k,0) -^E P 'V( £ P+k- £ p) 



-n(k) 



pV"p+k 



'P > (E p+k -E p ) 



+0(l/z 4 ) 

iSn(k, 0) iSn(k, 0) ^ay ^p n r 



Xo(k) 



■fi(k) 



+0(l/z 4 ). 



(41) 
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One can see that the coefficient of the term proportional 
to 1/z 2 vanishes, as is the case with Eq. (I39[) . 
As the initial condition, we set 

tfn(k,0) = C(k)F(k), (42) 

where, the function C(k) is given by 

C(k) 

duj JZo dt I d 3 x{[5n(x, t),Sn(0, 0)]) e ,e- ik V-* 



2tt 



(43) 



Substituting it into Eq. I|41(l . we find that the first term 
is identical with that of Eq. I|39(l . On the other hand, we 
can see the following relation by comparing the second 
terms: 



<n(x')) e9 = C(k) 



Xo(k) 



n(k). 



(44) 



When we apply the free gas approximation to the expec- 
tation value of Eq. I|43|) . the initial condition is given 
by 



C(k)=/?xo(k). 



(45) 



Substituting this expression into Eq. I|44|l . we obtain the 
number density at the equilibrium state, 



(0))e 9 -^E" P ^ k )- 



(46) 



The r.h.s. of the equation has the k dependence. It fol- 
lows that the number density converges to an inhomoge- 
neous distribution instead of a thermal equilibrium distri- 
bution. However, it should be noted that the breaking of 
the sum rule is small for low k because linu^o Q(k) = 1- 
Thus, this approach is still available for describing the 
dynamics associated with long distance scales [13 ■ 



V. MARKOV APPROXIMATION 

The derived diffusion equation is still an integrodif- 
ferential equation. To discuss the behavior of the coarse- 
grained equation, that is, causal or acausal, we derive the 
local equation by employing the Markov approximation. 

First of all, we should notice that the memory term 
can be separated into two terms (2jJ 0] , 

d f* 

—Sn(k,t) = - / dTn 2 Ck,t-T)Sn(k,T) 
at J 

- [ dT$(k,t-T)Sn(k,T), (47) 
Jo 

where the frequency function f2 2 (k, t) and the renormal- 
ized memory function ^(k, t) are defined by the imagi- 
nary part and the real part, respectively, 

fl 2 (k,t) = J ^zIm[r L (k,- l ^ + e)]e-^ t , (48) 



$(k,t) - J ^Re[r L (k, 



e)]e 



(49) 



It is worth notifying that the time derivative of the 
number density always vanishes at the initial time, 
d{5n{k,t))/dt\ t=0 = 0. 

From the final value theorem, we can see the tempo- 
ral behavior of the two functions; the frequency function 
converges to a finite value and the renormalized memory 
function vanishes at late time. Thus, we assume that the 
renormalized memory function relaxes rapidly and van- 
ishes, while the frequency function converges to finite val- 
ues depending on temperatures and chemical potentials 
after short time evolution. Actually, the memory func- 
tion of the chiral order parameter behaves as is assumed 
above 0] . Because we are interested in the slow dy- 
namics associated with the macroscopic time scale, we 
ignore such fast variations. Then, the frequency function 
is approximately given by a time-independent constant 
and the time dependence of the renormalized memory 
function is replaced by the Dirac delta function, 



Vt 2 (k,t) ps lim n 2 (k,t) = L> k k 2 , 



(50) 



$(k,i) ps 2S(t) / dr$(k,r) = — S(t). (51) 

JO T k 

As a result, the integrodifferential equation is approx- 
imately given by 



d 
dt 



f 4 1 

5n(k, t) = -L> k k 2 / drSn{k,T) Sn(k,t), 

Jo r k 



where 



2Q(k) f 2 _ 
m/W X0 (k) Jo PPnp> 



- = |k| Xo (k)n(k) 



1 



(2 



(52) 



(53) 



(54) 



Here, we introduced the momentum cutoff A. It should 
be noted that -D k is defined also by the final value theo- 
rem, 



£> k k 2 = limsP^k, s). 

s— >0 



(55) 



This definition gives completely the same expression as 
Eq. 

These expressions are further simplified in the low mo- 
mentum limit, 



£ k ps D = 



1 Jo dpp 2 n p 
n f3 f Q A dpp 2 rip (1 - rip) 
— pa 2 J^dpp 2 n p (l - n~) 
T k ~ "»7T J Q A dppn p (l-np) ' 



(56) 



(57) 
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For later convenience, we multiply rt for both sides. 
Finally, the integrodiffcrcntial equation is reduced to the 
following Markovian diffusion equation: 

rra<5n(k, t) + ^-Jn(k, t) + D c k a k 2 5n(k, t) = 0, (58) 
at z at 



where 



D c k a = Dr k . 



(59) 



Solving the differential equation, we should employ the 
constraint for the initial conditionl38l. 



j t Sn(k,t) 



= 0. 



(60) 



One can easily see that the equation is nothing but the 
causal diffusion equation. The diffusion coefficient £> k a 
corresponds to the diffusion constant of the usual dif- 
fusion equation, and represents the relaxation time. 
One can easily see that the Markovian diffusion equation 
is reduced into the form of the acausal diffusion equation 
by setting Tk = 0. 

The difference between Eq. (|2J and Eq. is the 

momentum dependence of the relaxation time. In Eq. 
(|2*|) , the number density decays with the same relaxation 
time even for the different momentum modes. On the 
other hand, our equation shows different relaxation time 
depending on the momentum. 



VI. DIFFUSION CONSTANT AND 
RELAXATION TIME 

If the coarse-grained dynamics is given by the acausal 
diffusion equation, the diffusion constant is expressed 
by the time correlation function of the number density 
(or the time correlation function of the number current 
density) |2lJ . In this section, we give the expressions of 
the diffusion coefficient and the relaxation time in terms 
of correlation functions. 

By using the initial condition given by the linear re- 
sponse of the external field F(k), the Laplace- Fourier 
transform of the causal diffusion equation (|58|l is pre- 
sented by 



5n LF {k,z)/F{k) 



(1 - izr^Cik) 



— Z Tk — IZ - 



(61) 



is related 



Note that the real part of 6n(k, z) / F(k) 
to the correlation function of the number density (see 
Appendix [BJl . 

Setting z = u) + it, we can find that the real part has 
the following relations; 



n(k,w) _ C(0)D 



^ok^o|k| 3 F{k) t 

, , „ n(k,u) £7(0) 
lim lim k Re-^-^ = 



(62) 
(63) 



where r = Umk->o(|k|Tk) and D is defined by Eq. 
H56JI . Thus, it is possible to express the diffusion co- 
efficient and the relaxation time in terms of the cor- 
relation function. However, the expressions are not so 
simple as the diffusion constant of the acausal diffu- 
sion equation. This is because our equation is based 
on the coarse-grained dynamics using the projection op- 
erator. In the POM, the projected microscopic vari- 
ables are the origin to cause the fluctuations and dis- 
sipations of variables associated with macroscopic time 
scales. Then, the diffusion coefficient is defined by the 
noise-noise correlation instead of the current-current cor- 
relation H IH HJ IH |H |H H3. As a mat- 
ter of fact, from Eq. (|L 1|) . the exact TC equation of the 
number density is given by 



d 
di 



Sn(x,t) — J dr J d 3 x'T(x — x' ,t — s)8n(x , s) 



+£(x,t). 



(64) 



Here, the last term represents the noise term that has 
been ignored until now and is defined by 



£(x,t) = Qe lLQt iLSn(x,0). 



(65) 



The relaxation function r(x — x', t — s) is expressed by 
the noise-noise correlation, 

T(x - x't - t') 

d 3 x"(£(x, t), £(x", t')) ■ (Sn(x", 0),6n(x>, 0))" 1 . 



(66) 

In the linear response theory, it is known that the dif- 
fusion constant is expressed by the correlation function of 
currents when the coarse-grained dynamics is assumed to 
be given by the acausal diffusion equation 21]. To derive 
the acausal diffusion equation after taking the Markov 
limit in the POM, the noise-noise correlation must be 
replaced by the current-current correlation, at least, in 
the low momentum limit. This is easily verified if we can 
approximately replace the coarse-grained time-evolution 
operator e li< 5* with the usual time-evolution operator 
e lLt (see Appendix lAl for details). As a matter of fact, 
this approximation means to drop \s s ) m the denom- 
inator of the memory function, and hence the memory 
function is approximately replaced by Eq. I|36l) . Then, 
one can easily show that the Markovian diffusion equa- 
tion is reduced to the acausal diffusion equation in the 
low momentum limit. It is normally assumed that the 
approximation is justified at least in the low momentum 
limit (2^1 ■ Actually, the famous conclusion that the diffu- 
sion constant is given by the current-current correlation is 
derived under the same assumption |39j| . However, xf(k) 
is not small even in the low momentum limit. This is the 
reason why we cannot obtain the acausal diffusion equa- 
tion in contrast to Ref . ■ This result further suggests 
that the current-current correlation does not necessarily 
give the definition of the diffusion constant. 
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VII. SUMMARY AND CONCLUDING 
REMARKS 

We applied the projection operator method (POM) to 
the non-relativistic model and derived the coarse-grained 
equation for the number density. The derived equation is 
an integrodifferential equation and contains the memory 
effect. In our model, the number density is a conserved 
quantity and there exists the sum rule associated with 
it. The usual acausal diffusion equation breaks the sum 
rule. On the other hand, the integrodifferential equation 
satisfies the sum rule in the low momentum limit. 

Second, we assumed that there exists the clear sepa- 
ration between microscopic and macroscopic scales, and 
employed the Markov approximation. The Markovian 
diffusion equation is characterized by the diffusion coef- 
ficient and the relaxation time as is the case with the 
causal diffusion equation. Thus, we can conclude that 
the causal diffusion equation can be derived by using the 
POM. However, it should be noted that the relaxation 
time of our equation depends on the momentum. 

To derive the acausal diffusion equation in the POM, 
as is discussed in Ref. |23 , we should approximately re- 
place the coarse-grained time-evolution operator with the 
normal time-evolution operator. Usually, it is assumed 
that the approximation is justified in the low momentum 
limit |22j|. However, in our calculation, we cannot em- 
ploy the approximation even in the low momentum limit. 
This result further suggests that the current-current cor- 
relation does not necessarily give the definition of the 
diffusion constant. 

These results are not particular to the model used in 
this paper. The causal diffusion equation is obtained 
also by applying the POM to the Nambu-Jona-Lasinio 
(NJL) model, which is the low energy effective model of 
quantum chromodynamics (QCD). This result may be 
important to discuss the QCD critical dynamics. Usu- 
ally, we simply assume acausal diffusion equations as 
coarse-grained equations of conserved variables associ- 
ated with macroscopic time and length scales near the 
critical points 0, 

El E E3- However, the coarse- 
grained equation does not necessarily have an acausal 
form [231 l24l lil ] . When we apply causal diffusion equa- 
tions instead of acausal ones, the QCD critical dynamics 
might be changed. 

In this calculations, we employed the perturbative ap- 
proximation. Thus, there may exist criticism that if we 
calculate without employing approximations, we may be 
possible to obtain acausal diffusion equations instead of 
causal ones. However, it is impossible because the exact 
calculation must satisfy the sum rule. 

Ichiyanagi also discussed the causal diffusion equation 
(more generally, the extended irreversible thermodynam- 
ics) based on the projection operator method 0, liq . 
However, the situation discussed by him is different from 
ours. He introduced a time-smoothed density matrix to 
define the projection operator. The time-smoothed den- 
sity matrix represents a non-equilibrium state and hence, 



the Kubo's canonical correlation of Eq. i|24|) is replaced 
by a non-equilibrium expectation value. In our calcula- 
tion, we consider the situation where the deviation from 
the equilibrium state is not so strong and the degrees 
of freedom associated with microscopic time scale has 
already reached the thermal equilibrium. This is the 
reason why we employ the Kubo's canonical correlation 
as the definition of the scalar product contained in the 
projection operator. On the other hand, when we em- 
ploy the time-smoothed density matrix as the definition 
of the scalar product, the microscopic degrees of free- 
dom stays in the non-equilibrium state and hence the 
derived coarse-grained equation describes the evolution 
of the fluctuations from the non-equilibrium state. 

Recently, the generalized projection operator method 
(GPOM) is developed to describe the pulse parameter 
dynamics |47^ . Although projection operators are intro- 
duced also in the GPOM, the GPOM is different from our 
method. The Mori projection operator satisfies the con- 
dition 10, and this plays an important role in deriving 
the TC equation (|I6|) . However, the projection operator 
introduced in the GPOM does not satisfy this condition 
because the phase degree of freedom is introduced. Thus, 
it is difficult to reproduce the GPOM in our projection 
operator method. 
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APPENDIX A: CONVENTIONAL DERIVATION 
OF DIFFUSION EQUATION IN THE 
PROJECTION OPERATOR METHOD 

In this section, we review the conventional derivation 
of a coarse-grained equation of the number density in the 
POM 22]. 

Substituting the Mori projection operator defined by 
Eq. (|2*2*|) into Eq. iJTT]). the exact TC equation is given 
by 

= J dr y"d 3 x'd 3 x"(e iQLT (5n(x),(5n(x')) 
■(5n(x'), (5n(x")) _1 • 5n{x", t - r) + £(x, t), 

(Al) 

The streaming term vanishes in this definition of the pro- 
jection operator. We assume that the memory function 
can be approximated as follows: 

(e iQLT Sn( x ),8h( x ')) = (e iQLT V x -SJ(x),V x >-5J(x')) 
« -vt(e iiT 5J(x),<5J(x')), (A2) 
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where <$J(x) is the current of the number density. In the 
first line, we used the equation of continuity, 



dt 



5n(x,t) + W -63(x,t) = 0. 



(A3) 



In the second line, the coarse-grained time-evolution op- 
erator e l ® Lt is approximately replaced by the usual time- 
evolution operator e lLt . 

After the Fourier transformation, the TC equation is 
given by 



lit 



Sn(k,t) = k 2 f dTr.jj(k,T)5n{k,t -t) 
Jo 

+e(M), (A4) 



where Tjj(k,t) is the Fourier transform of the memory 
function, 



cfx'(e' Lt <5J(x),,5J(x')) • (5n(x'),Sn(x")Y 



(AS) 



We assume that Tjj(0,t) is finite after taking the 
Markov approximation. Finally, the coarse-grained equa- 
tion at low k is given by 

^-Sn(k, t) = k 2 D aca Sn(k, t) + £(k, t), (A6) 
dt 

where the diffusion constant is defined by 



POO 

D aca = / drYjj{0,T). 
Jo 



(A7) 



This is the acausal diffusion equation, and hence it is 
often claimed that acausal diffusion equation can be de- 
rived in the POM. However, as we have seen, to derive 
the acausal diffusion equation, we should apply the ap- 
proximation (|A2|) and it is not applicable to the non- 
relativistic Hamiltonian used in this paper. 



APPENDIX B: CORRELATION FUNCTIONS OF 
NUMBER DENSITY 

First of all, we summarize the general properties of the 
correlation functions of the number density following the 
discussion given by Kadanoff and Martin 15]. This dis- 
cussion is of assistance when we investigate the validity 
of the coarse-grained equation obtained by applying the 
POM. 

In the linear response theory, the expectation value of 
an arbitrary operator (0(t)) that is dynamically induced 
by the external field is given by 



(0{x,t)) - (0(x,t)) eq =i I ds([H ex (s),0(x,t)}) eq , 

(Bl) 



where ( ) eq means to take a thermal expectation value 
with a Hamiltonian H, and 0(x,t) = e iHt O(x)e~ iHt . 

We are interested in the number density induced by 
the external field -F(x, t). For this purpose, we substitute 
O(x) = n(x) and H ex {t) = — J d 3 xn(x,t)F(x,t) into Eq. 
ill. Then, we obtain 



(Sn(x, t))=i dt'([n(x, t), n(x', t')]) eq F{x' , i').(B2) 



Here, the external field has the following time depen- 
dence: 



F(x,t) 



F(x)e et t < 



t > 0. 

Now, we introduce a function C"(k, uS) as follows; 

)eq 



(B3) 



<[n(x,*),n(x',t')]>e 



du> f dk 



2tt J (2tt) 3 



C"(k,c t ;)e ik ( x - x ')e-^(*- t '). 



(B4) 



The C"(k, lu) is real and an odd function of the frequency 

LJ. 

Substituting Eq. l|B4p into Eq. I|B2|). we obtain 

(6n(x, t)} 



dtu f d 3 k ,C"(k, 



U3) 



2tt J (2tt) 3 



F(k) " e tkx e- iui (t > 0). 



(B5) 



We further define the Laplace-Fourier transform of the 
number density, 



&n L *%z) = d A x dte lzt {Sn(k,t)) 



du C"(k,u) 
2ni uj(uj — z) 



F(k). 



Here, we substitute Eq. (|B5|I . 

The dynamic susceptibility is defined by 



(B6) 



C(k, z) 



duj' C"(k,u/) 
2ir uj' — z 



(B7) 



Setting z = to + ie, the dynamic susceptibility is decom- 
posed into the real part and the imaginary part; 



C(k, w + ie) = C"(k, u) + ^C"(k, uj), (B8) 



where 



C'(k,uj)=P J 



dtu' C"(k,w') 

27T U>' — UJ 



(B9) 



This is a Kramers-Kronig relation. Because the function 
C"(k, ui) is an odd function of uj, we obtain the following 
relation from the Kramers-Kronig relation, 



C(k) = C(k,uj = 0) = J 



duj' C"(k,uj') 
~2ir uj' 



(BIO) 
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We can derive a sum rule with the help of the equation 
of continuity, 



0_ 

Of 



n(x,t) +V- J(x,t) =0. 



(Bll) 



Operating the time derivative to Eq. i|B4|) and applying 
the equation of continuity, we obtain 



d_ 

Of 



<Hx,i),n(x',i')]> 



eq 



= -([V-J(x,t),n(x',t')]) 



eq 



dw f dk 



2tt J (2tt) : 



wC*(k,u)e 



£k(x— x') — io;(t— i') 



(B12) 



In the Schrodinger field discussed in this paper, the 
current operator J(x, t) is given by 

J(x,t) = ^-(V^(x,i) • V(x,i) -^( Xji )vV»(x,t)). 

(B13) 

Then, one can easily calculate the equal time commutator 
of the current and the number density, 

<[J(x, t), n(x', t)]) eq = --V<5(x - x')(n(x')} eg . (B14) 



This equation implies 

<iu> /" dk 



2tt 7 (2?r) : 



^C"(k,w)e tk ( x - X ') 



1 



= -V 2 <5(x-x')(n(x' = 0)) e „ 
and hence we can obtain the following f-sum rule, 



(B15) 



2^ 



1 



wC"(k,w) = — k 2 (n(0)} 



(B16) 



By expanding Eq. (|B6|I for large values of z and using 
the relation and the sum rule, we obtain 

5n LF (k,z)/F(k) = -C(k) + 1 / ^C"(k,u,) 
z z J 2n 



* „/. s i 1 



3 , _W<7"(k, W ) + 



= -C(k) + 3 -k 2 (n(0)) e(? + 0(l/z 4 ). 

(B17) 



Z J TO 



1. Acausal diffusion equation 

Here, we assume that the time-evolution of the num- 
ber density follows the acausal diffusion equation. Then, 
from the Fick's law, the current is given by 



J(x,i) = £)Vn(x,t), 



(B18) 



where D is the diffusion constant. Substituting it into 
the equation of continuity, we obtain the acausal diffusion 
equation, 

d 

— n(x, t) = L>V 2 n(x, t). (B19) 

After the Laplace-Fourier transformation, the number 
density is given by 

n(k,0) 



n LF (k,z) 



Dk 2 ' 



(B20) 



where, n(k, 0) represents an initial value of the number 
density. When the initial value is induced by the external 
field F(x, t) defined in the preceding section, we can set 



i(k,0) = C(k)F(k). 



(B21) 



Expanding the Laplace- Fourier transform for large values 
of z, we have 



n iF (k,z) CCk) C(k)Dk 2 
= i t. 



F(k) 



(B22) 



This expression has the term proportional to l/z 2 , while 
the corresponding term vanishes in Eq. (|B17|) . Further- 
more, C"(k,a->) in the diffusion equation is 

C(k)Dk 2 oj 



C"(k,u) 



uj 2 + {Dk 2 f 



(B23) 



This expression fails to satisfy the sum rule (|B16|) . Thus, 
if the coarse-grained dynamics is approximated by the 
acausal diffusion equation, the time-evolution completely 
breaks the sum rule. 



2. Causal diffusion equation 

Kadanoff and Martin pointed out that we should intro- 
duce a relaxation time to satisfy the sum rule[l5|. Then, 
the Fick's law is modified as follows; 



d T / \ 



— J(x, t) - DVn(x,f). 



(B24) 



Substituting this equation into the equation of continuity, 
we obtain 

dt 2 



ld_ 

rdt 



n(x,t) = 0. 



(B25) 



After using the initial conditions 

n(k,0)=C(k)F(k), 

-0, 



m n(x,t) 



(B26) 
(B27) 



t=o 



the Laplace-Fourier transform of the number density is 
given by 

n LF (k,z) C(k){l-izT) 



F(k) -iz + Dk 2 ~ tz 2 

.C(k) Dk 2 C(k) 



(B28) 
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In this expression, we do not have the term proportional 
to l/z 2 and the equilibrium value of the number den- 
sity is given by (n(0)) eq — mC(k)D/T. The equilibrium 
number density should not depend on the momentum and 
hence the expression of (n(0)) eq looks inconsistent. How- 



ever, the diffusion equation is the coarse-grained equation 
and it is valid only in the low momentum limit. Then, the 
function C(k) is approximately given by C(0). Thus, one 
can conclude that the causal diffusion equation is consis- 
tent with the sum rule in the low momentum limit . 



[1] C. Cattaneo, Atti. Semin. Mat. Fis. Univ. Modena 3, 83 
(1948). 

[2] I. Miiller, Z. Phys. 198, 329 (1967). 
[3] W. Israel, Ann. Phys. (N.Y.) 100, 310 (1976). 
[4] W. Israel and J. M. Stewart, Ann. Phys. (N.Y.) 118, 341 
(1979). 

[5] I. S. Liu, I. Miiller, and T. Ruggeri, Ann. Phys. (N.Y.) 

169, 191 (1986). 
[6] R. Geroch and L. Lindblom, Ann. Phys. (N.Y.) 207, 394 

(1991). 

[7] P. Kostadt and M. Liu, Phys. Rev. D 62, 023003 (2000). 
[8] D. Jou, J. Casas- Vazquez, and G. Lebon, 

Rep. Prog. Phys. 51, 1105 (1988). 
[9] D. Jou, J. Casas- Vazquez, and G. Lebon, 
Rep. Prog. Phys. 62, 1035 (1999). 
[10] A. Muronga, Phys. Rev. Lett. 88, 062302 (2002); 89, 

159901(E) (2002). 
[11] A. Muronga, Phys. Rev. C 69, 034903 (2004). 
[12] A. Muronga and D. H. Rischke, nucl-th/0407114 
[13] M. Abdel-Aziz and S. Gavin, Phys. Rev. C 70, 034905 
(2004). 

[14] W. L. Kath, Physica D 12, 375 (1984). 
[15] L. P. Kadanoff and P. C. Martin, Ann. Phys. 24, 419 
(1963). 

[16] J. C. Maxwell, Philos. Trans. R. Soc. London 157, 49 
(1867). 

[17] S. Machlup and L. Onsager, Phys. Rev. 91, 1512, (1953). 

[18] H. Grad, in Hundbuch der Physik, edited by S. Flugge, 
(Springer, Berlin, 1958) Vol. 12. 

[19] The problem of causality can be solved also by introduc- 
ing nonlinear terms fl4|l . 

[20] As a matter of fact, the deviation from the acausal dif- 
fusion equation has been discussed in, for instance, ther- 
mal diffusion processes, spin diffusion processes and so 
on 0. 

[21] R. Kubo, M. Toda, and N. Hashitsume, Statistical 
Physics II ( Springer- Verlag, Berlin, 1983). 

[22] E. Fick and G. Sauermann, The Quantum Statistics of 
Dynamic Process (Springer- Verlag, Berlin, 1983). 

[23] T. Koide and M. Maruyama, |nucl -th/0308025l 

[24] T. Koide and M. Maruyama, JNucl. Phys. A 742, 95 
(2004). 



[25] S. Nakajima, Prog. Theor. Phys. 20, 948 (1958). 

[26] R. Zwanzig, J. Chem. Phys. 33, 1338 (1960). 

[27] H. Mori, Prog. Theor. Phys. 33, 423 (1965). 

[28] N. Hashitsume, F. Shibata and M. Shingu, 

J. Stat. Phys. 17, 155 (1977). 
[29] F. Shibata, Y. Takahashi and N. Hashitsume, 

J. Stat. Phys. 17, 171 (1977). 
[30] F. Shibata and T. Arimitsu, J. Phys. Soc. Jpn 49, 891 

(1980). 

[31] C. Uchiyama and F. Shibata, Phys. Rev. E 60, 2636 
(1999). 

[32] T. Koide, M. Maruyama and F. Takagi, 

Prog. Theor. Phys. 101, 373 (1999). 
[33] T. Koide and M. Maruyama, Prog. Theor. Phys. 104, 

575 (2000). 

[34] T. Koide, Prog. Theor. Phys. 107, 525 (2002). 

[35] J. Rau and B. Miiller, Phys. Rep. 272, 1 (1996). 

[36] U. Balucanai, M. Howard Lee, and V. Tognetti, 

Phys. Rep. 373, 409 (2003). 
[37] As a matter of fact, even for the ori gina l causal diffusion 

equation [llSliUHSEim El 111 , we need to take 

the low momentum limit to preserve the sum rule [l5|. 

See Appendix iBl 
[38] The constraint derived from Eq. 1521 is 

d(«5n(k, t))/dt\ t=0 — — Sn(k, 0)/rk. However, as is 

discussed below Eq. 1491 . we should employ Eq. (1601 . 

This difference is caused by the Markov approximation. 
[39] L. E. Reichl, A Modern Course in Statistical Physics, 

(University of Texas Press, Austin, TX, 1980). 
[40] K. Rajagopal and F. Wilczek, Nucl. Phys. B399, 395 

(1993). 

[41] D. T. Son and M. A. Stephanov, hep-ph/0401052 

[42] K. Ohnishi, K. Fukushima anT K~ Ohta, 

|micT- th/04090l6l 

[43] T. Koide, hep-ph/0411207 

[44] T. Koide, G. Krein and R. O. Ramos, in preparation. 
[45] M. Ichiyanagi, J. Phys. Soc. Jpn, 59, 1970 (1990). 
[46] M. Ichiyanagi, Prog. Theor. Phys. 84, 810 (1990). 
[47] K. Nakkeeran and P. K. A. Wai, Opt. Commun. 244, 377 
(2005). 



